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We consider (do + 2)-dimensional configurations with global strings in the two extra dimensions 
and a flat metric in do dimensions, endowed with a warp factor e 27 depending on the distance I 
from the string center. All possible regular solutions to the field equations are classified by the 
behavior of the warp factor and the extra-dimensional circular radius r(l). Solutions with r — ► oo 
and r — * const > as I —* oo are interpreted in terms of thick brane world models. Solutions 
with r — > as I — > l c > 0, i.e., those with a second center, are interpreted as either multi-brane 
systems (which is appropriate for large enough distances l c between the centers) or as Kaluza- 
Klein type configurations with extra dimensions invisible due to their smallness. For the case of 
the Mexican-hat symmetry breaking potential, we build a full map of regular solutions on the 
(e, r) parameter plane, where e acts as an effective cosmological constant while T characterizes 
the gravitational field strength. The trapping properties of candidate brane worlds for test scalar 
fields are discussed. Good trapping properties for massive fields are found for models with growing 
warp factors. Kaluza-Klein type models are shown to possess nontrivial warp factor behaviors, 
leading to matter particle mass spectra which seem promising from the standpoint of hierarchy 
problems. 

PACS: 04.50.+h, 11.27.+d 
1. Introduction 



The multidimensional gravity concept, tracing back to the pioneering papers by Kaluza and Klein [1], 
initially assumed that the extra dimensions remain unobservable due to their extreme smallness. 
Another class of multidimensional theories has been put forward in the 80s, on the basis of the 
idea that we live on a distinguished surface (brane) embedded in a higher- dimensional manifold, 
called the bulk [2]. This idea has recently become very popular in attempts to find an approach 
to a number of fundamental physical problems. The brane-world concept is broadly discussed in 
connection with the recent developments in supersymmetric string/M-theories [3]. The simple RSI 
model [4] has continued numerous attempts [5] to find the origin of the enormous hierarchy of 
energy/mass scales observed in nature, which is a long-standing problem in particle physics. In 
astrophysics and cosmology, there are attempts to explain the dark matter and dark energy effects, 
to describe black holes and possible variation of fundamental constants, the CMB anisotropy etc. 
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A great variety of brane world models may be found in the literature: branes in 5 or more 
dimensions, single or multiple branes, flat or curved branes, flat or curved bulk, compact or non- 
compact bulk (i.e., large or infinite extra dimensions), thin or thick branes, various symmetries of 
both bulks and branes, different kinds of matter forming the brane etc. 

In our view, the most natural physical idea leading to the emergence of distinguished surfaces 
in the space-time manifold is the idea of a phase transition with spontaneous symmetry breaking 
(SSB) which has already led to great success in many areas of physics and cosmology. In other 
words, it is reasonable to consider the brane world as a result of a phase transition at a very early 
stage of the Universe evolution. The existing macroscopic theory of phase transitions with SSB 
allows one to consider the brane world concept self-consistently and maximally escape the influence 
of model assumptions even without a detailed knowledge of the nature of physical vacuum. A 
necessary consequence of such phase transitions is the appearance of topological defects. 

Let us recall that a decisive step toward cosmological applications of the SSB concept was made 
in 1972 by Kirzhnits [6]. He assumed that, as in the case of solid substances, a symmetry of a field 
system, existing at sufficiently high temperatures, could be spontaneously broken as the temperature 
falls down. The first quantitative analysis of the cosmological consequences of SSB was given by 
Zel'dovich, Kobzarev, and Okun' [7]. 

The properties of global topological defects are generally described with the aid of a multiplet 
of scalar fields playing the role of an order parameter. If a defect is to be interpreted as the origin 
of a brane world, its structure is determined by the self-gravity of a scalar field system and may be 
described by a set of Einstein and scalar equations. 

This approach to the brane- world concept has been employed by many authors for construction of 
thick branes in five (see, e.g., [8-14]) and more (see, e.g., [15-18] and references therein) dimensions. 

In our previous papers [13,14,17,18], we have analyzed the gravitational properties of candidate 
(thick) brane worlds with the do -dimensional Minkowski metric and global topological defects in 
d\ + 1 extra dimensions. Our general formulation covered such particular brane (domain 

wall) in five-dimensional spacetime (one extra dimension), a global cosmic string with winding 
number n = 1 (two extra dimensions), and global monopoles (three or more extra dimensions). We 
restricted ourselves to Minkowski branes since most of the existing problems are clearly seen even in 
these comparatively simple systems; on the other hand, in the majority of physical situations, the 
intrinsic curvature of the brane itself is much smaller than the curvature related to brane formation, 
and therefore the main qualitative features of Minkowski branes should survive in realistic curved 
branes. 

Our treatment differed from many others (e.g., [15,16]) in that we have considered all kinds of 
regular solutions to the corresponding field equations, including those with growing warp factors, 
and have emphasized good trapping properties of the latter. 

We have shown, in particular, that there are as many as seven classes of regular solutions of 
the field equations describing global strings and monopoles in extra dimensions; two of them exist 
for monopoles only while the other five are found for both strings and monopoles. Some of these 
configurations have exponentially growing warp factors ( e 27 in the metric ([I]), see below) at large 
distances from the core. They are shown to trap linear test scalar fields of any mass and momentum. 
Others, ending with a flat metric, have a warp factor tending to a constant value, determined by 
the shape of the symmetry breaking potential. They are also shown to trap test scalar fields with 
masses bounded above by a value depending on the particular parameters of the topological defect. 

Though the general classification of [17, 18] covers all possible regular configuration, the impor- 
tant question of location of different solutions in the space of physical parameters remained open. 
One of the goals of this paper is to answer this question for the particular case of a global string in 
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two extra dimensions and a Mexican-hat symmetry breaking potential. The problem then contains 
two essential physical parameters: e, the dimensionless cosmo logical constant, and T , characterizing 
the gravitational field strength. The border lines in the (e, T) plane, separating different classes of 
regular solutions (those extending to infinite circular radii r, those with a cylindrical geometry far 
from the center and those with with two regular centers), are found numerically, and the asymptotic 
dependences e(T) as T — > and r — > oo are derived analytically. 

Another goal is to give a more complete description of the configurations of interest described by 
these solutions. Thus, we describe the trapping properties of different thick brane world models for 
classical particles, scalar fields and gravity. We also argue that one of the classes of regular configu- 
rations, those with two centers, can lead to promising models with nontrivial particle mass spectra. 
The point is that the warp factors of such configurations can have several minima at different levels, 
where test particles and fields may be gravitationally trapped with different energies. Though, in 
this case, one should abandon the brane world concept and admit that the extra dimensions are 
invisible due to their smallness, i.e., interpret the solutions in the spirit of Kaluza-Klein theories. 

There is a growing number of publications devoted to brane worlds with two extra dimensions, 
see [19] and references therein. In many cases the results are obtained numerically using simplified 
models with specially chosen sets of parameters. In our macroscopic approach, based on the theory 
of phase transitions with SSB, we try to reduce the influence of model assumptions to a minimum 
and to cover the full range of opportunities. Our main result, the full map of regular solutions for a 
system with the Mexican-hat symmetry-breaking potential V(4>) , should probably retain its basic 
qualitative features for other potentials with a similar disposition of extremum points. 

The paper is organized as follows. Sec. 2 outlines the problem setting, including the geometry, 
field equations and boundary conditions providing space-time regularity. Sec. 3 describes the sim- 
plest solutions with a constant scalar field, needed for comparison in what follows. Sec. 4 is central 
in the paper: we give a general description and classification of possible regular solutions on the 
basis of our previous work [17,18] and present a map showing the location of different solutions with 
the Mexican hat potential in the parameter space of the problem. In Sections 5 and 6 we describe 
some further details of the solutions and give an analytic derivation of the asymptotic behavior of 
the curves drawn in the map. Sec. 8 discusses the trapping properties of thick branes described by 
the above solutions. Sec. 9 outlines the properties of configurations with two centers, and Sec. 10 is 
a conclusion. 



2. Problem setting 

2.1. Geometry and regularity conditions 

Our main interest here is 6D space-time with a cosmic string in the two extra dimensions. Let us, 
however, begin with a more general geometry: consider a (D = d + d\ + 1) -dimensional space-time 
with the structure M. d ° x R u x E dl and the metric 

ds 2 = e 2 ^ u \^dx»dx v - ( e 2a ^du 2 + e 2l3(u) dfl 2 ) . (1) 

where r]^ = diag(l, — 1, — 1) is the do -dimensional Minkowski metric (do > 1), dQ is a linear 
element on a d\ -dimensional unit sphere E> dl , a, (3 and 7 are functions of the radial coordinate 
u with the definition domain M u C 1 to be specified later. We will also use the notation r = e 13 
where r is the spherical (circular for d± — 1) radius. 

The Riemann tensor R AB cd is diagonal with respect to pairs of indices and has the nonzero 
components 
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R ab cd = (e-^-e- 2a (3' 2 )5 ab cd , 
R ua ub = -5 a b ^- a {^- a p')\ 

n = -«e" 2 ^' (2) 

where = 5^5% — 5^5 u p and similarly for other indices. The indices /i, v, ... correspond to d - 
dimensional (physical) space-time, a,b,... to d\ angular coordinates on § dl , and A,B,... to all D 
coordinates. 

A necessary condition of regularity is finiteness of all algebraic invariants of the Riemann tensor. 
In our case it is sufficient to deal with the Kretschmann scalar K = R a ^dR CL ab si nce it is a sum 
of squares of all nonzero R AB cd • Hence, in regular configurations all components of the Riemann 
tensor (jSJ) are finite. 

In the Gaussian gauge a = 0, u = I being the proper distance along the radial direction, the 
regularity conditions at r > look very simple: 

(3' , f3" , 7', 7" should be finite. (3) 

The regularity conditions at the center r = follow from finiteness of the Riemann tensor compo- 
nents R ab cd and coincide with the regular center conditions in usual static, spherically symmetric 
metrics. In terms of an arbitrary u coordinate, the regular center conditions are 

7 = 7c + 0(r 2 ) e p - a \(3'\ = l + 0(r 2 ) as r^O. (4) 

The latter condition means a correct (= 2ir) circumference to radius ratio, or, equivalently, dr 2 = 
dl 2 ; 7 C is a constant which can be set to zero by rescaling the coordinates . 

The string case d\ = 1 has a specific feature: there is only one angular coordinate, therefore 
S ab cd = 0, hence R ab c d = 0. However, a conical singularity is possible (i.e., an angular deficit, 
dr 2 < dl 2 , or excess, dr 2 > dl 2 ), which is a pointwise, delta-like curvature peak over this zero level, 
as in the case of an ordinary cone top. Its existence actually means that there is some pointlike 
object with respect to the two extra dimensions, or a thin brane in the space-time as a whole. 

Below we will consider entirely regular configurations, excluding, among others, conical singu- 
larities. 



2.2. Topological defects. Field equations 

A global defect with a nonzero topological charge can be constructed with a multiplet of d± + 1 
real scalar fields (p k , in the same way as, e.g., in [17]. It comprises a "hedgehog" configuration in 

R u x § dl : 

(f) k = <f)(u)n k {x a ), 

n k is a unit vector in the d\ + 1 -dimensional Euclidean target space of the scalar fields: n k n k = 1. 
The total Lagrangian of the system is taken in the form 

L = v~2 + lg AB dA<f> k d B <f> k - V(<f>), (5) 

where R is the .D -dimensional scalar curvature, x 2 is the D -dimensional gravitational constant, 
and V is a symmetry-breaking potential depending on <p 2 (u) = a a . 

The case d\ = , with only one extra dimension, is a flat domain wall. Regular thick Minkowski 
branes supported by scalar fields with arbitrary potentials were analyzed in [13,14]. 
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The case d± = 1 is a global cosmic string with the winding number n = 1, to be discussed 
here in detail. In case d\ > 2 we have a global monopole in the extra space-like dimensions, see, 
e.g., [15-18,20,21]. 

Let us write down the scalar field equation and three components of the Einstein equations for 
such systems in the Gaussian gauge a = 0, u = I (the prime denotes d/dl): 

<P" + (d i + - d 1 e" 2/3 

7 " + d ol ' 2 + d 1 p' 7 > 

f3" + d p'j' + d 1 f3' 2 
(dt/3' + d ol ') 2 - d ol ' 2 - d^' 2 
Any three of the above four equations are independent, and the fourth one is their consequence. 



dV 



D-2 ' 



(d 



l_A;V)e- 2/3 



2x 2 



-V, 



D-2 

x 2 (0' 2 - 2V) + d x e-^idx - 1 - x 2 2 ) 



(6) 

(7) 

(8) 
(9) 



2.3. Boundary conditions and fine-tuning relations 

The metric can be rewritten in the form 

ds 2 = e 2 ^ l) r]^dx t "dx u - dl 2 - r 2 {l)dQ 2 , (10) 

where r = e 13 is the spherical radius. Assuming that there is a regular center (r = 0), we put there 
/ = without loss of generality, and we will classify the configurations under study by the limiting 
value of r(l) (infinite, finite, or zero) at the largest or infinite values of /. 

In the general monopole case, the regular center requirement leads to the following boundary 
conditions for Eqs. ©-([9]) at / = : 

0(0) = 7 '(0) = r(0) = 0, r'(0) = 1. (11) 

The system (J6j) — (J9J) does not contain 7 but only its derivatives. It is convenient for numerical 
integration to work with Eqs. (El) — (EI) solved with respect to the second-order derivatives and consider 
([9]) as their first integral. 

We thus have four boundary conditions (|TT|) for the (effectively) fourth-order set of equations. It 
might seem that we must obtain a unique solution. However, this is not the case since / = , being 
a singular point of the spherical coordinate system (not to be confused with a space-time curvature 
singularity), is also a singular point of our set of equations. As a result, our set of equations admits 
an additional freedom of choosing <f)'(0) ; or, instead, we may use the requirement of global regularity 
to obtain a unique solution. 



Infinite extra dimensions 

If the solution is defined in the interval < / < 00, then the lacking boundary condition can be 
taken as — > const as / — > 00 , or 

0'(oo) = 0. (12) 

In general, when the scalar field starts from a maximum of the potential and ends at a finite value 
of 0, the five boundary conditions ( TTTT) and fll2jl uniquely determine a nontrivial solution to our 
field equations. Its existence determines a certain area in the space of parameters of the problem 
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without a priori fine tuning. An asymptotic analysis at I — > oo shows that in this general case f3(l) 
is a linearly growing function as I — > oo, and r'(oo) > 0. 

r'(oo) = is the special case when r tends to a finite constant at I — > oo. The solution 
then terminates at a slope of the potential rather than its minimum. The supplementary condition 
r'(oo) = seems to require a fine tuning relation between the free parameters of the problem. But 
we shall see that it is not quite so. An analysis of solutions with = O — const shows that there is 
also an area in the parameter space where the condition (|12p is fulfilled automatically, and solutions 
with r'{l) — > as I —* oo exist without any fine tuning. 

Two centers 

It can happen that the integral curves of a solution terminate at some finite value l c with r(l c ) = 
and 0(/ c ) = 0, which is one more center. Of interest for us are configurations in which this second 
center I = l c is also regular. Then the same four conditions (fTTI) should hold at I — l c too. Two of 
them can be fulfilled by choosing the values of <f>'(0) and l c . Two others can only be fulfilled by a 
proper choice of free (input) parameters of the problem, e.g., those of the potential (if any) and the 
cosmo logical constant. 

In the special case of symmetry between the center^], the input parameters are connected by 
only one fine-tuning relation. In this case, the boundary conditions at the second center are satisfied 
automatically, and the existence of a regular solution is provided by three conditions of smoothness 
at the middle (equator) point Z cq = l c /2: 

0'(W = 7'(W> = r'(l eq ) = 0. 

Two of these three conditions determine the values of l c and <p'{0), and the remaining one requires 
fine tuning of the input parameters of the problem. 

The fine tuning could be avoided at the expense of admitting conical singularities (in the string 
case) at the two centers. For symmetric solutions the three smoothness conditions at the equator 
could be satisfied by choosing (f)'(0), r'(0) and l c . In the general case of nonsymmetric centers, the 
four conditions of the second center could be satisfied by appropriately choosing 4>'(0), r'(0), r'(l c ) 
and l c . 



1 Eqs. are invariant under translations I — ► I + Iq and under reflections Iq + I — > Iq — I, <fi — > — (f>. This 

invariance gives rise to existence of regular solutions with two centers, which are symmetric against the middle 
point [17]. 

Furthermore, a solution with two regular centers defined in the interval (0, l c ) can be symmetrically extended to 
the next interval (l c ,2l c ) and further on, thus leading to a periodic solution for I £ R [22]. The metric remains 
regular everywhere; however, the points of contacts (0, ±l c , ±2i c , ...) are geometrically ambiguous: each of them 
belongs to two adjacent manifolds simultaneously. 

If one still believes in the reality of such systems, one can note that the spectrum of a low-energy particle in a 
perfectly periodic potential has a conductivity zone, allowing free propagation and thus making the extra dimension 
observable in principle. However, if the conductivity zone is very narrow, then even small perturbations should lead 
to localization of a particle. This interesting new possibility is worth a special consideration. 
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2.4. String equations 

In the string case d% — 1 to be considered here in detail, the field equations take the form 

0" + (cio7 / + /5 ,2 ¥-e- 2/3 0=|^, (13) 

T " + d 0l ' 2 + = -^-V, (14) 

{3" + d /3V + f3' 2 = -(x 2 2 ) - *£v, (15) 
2d 7'/3' + d (d - l)i 2 = x 2 (0' 2 - 2V) - *?<t> 2 e" 2/3 . (16) 
For numerical examples, we will use the so-called Mexican hat potential 
1 r / ch 2 \ 2 i 

V = -\oV 4 e+(l-^) . (17) 



4 



r/ 2 



The parameter e plays the role of a cosmological constant added to the conventional Mexican 
hat potential (the "hat" is thus moved up or down). To pass to dimensionless quantities, we put 
A ?7 2 = 1, so that lengths will be measured in units of l/(\/Ao?7) which in many cases has the 
meaning of the string core radius^). 
The potential then takes the form 

V=\v 2 [e + (l-m /:=*• (18) 

The remaining free parameters that control the system behavior are 

d , e, and T := x 2 ^ 2 . (19) 

We will use d = 4 in computations. The parameter T characterizes the gravitational field strength. 

In the further description of different regular solutions to our field equations, we begin with 
the simplest solutions in which = 0* = const. They are not string solutions but are helpful for 
comparison. 



3. Solutions with = const 

If = 0* = const , we are actually dealing with vacuum field equations for the metric (jTU!) with the 
cosmological constant A = x 2 V^(0*). 

The scalar field equation (1131) reduces to 

0*/r 2 = -W*), V+ := dV/d<j). (20) 

This leads to two kinds of solutions: one exists in the case V^(0*) = 0, 0* = 0, the other corresponds 
to ^(0*) 7^ (i.e., is "frozen" on a slope of the potential), and in this case we should put 
r = const. 

For the potential (jTTj) . Eq. p)]) gives 

/*(r- 2 -l + / 2 ) = 0, (21) 

so that either /* = 0*/r] = or / 2 = 1 — 1/r 2 . 

2 From the very beginning, we put c = H = 1 , so that all quantities are measured in appropriate powers of length 
[£} . Then some relevant dimensionalities are [V] = [t~ D ] , [<j> 2 } = [rj 2 ] = \x~ 2 \ = [^ 2_I? ] . 
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3.1. Solutions with = 

The trivial regular solutions with the order parameter = describe configurations with a higher 
symmetry, which can become spontaneously broken into a structure with a topological defect. 
In this case the metric obeys the equations 

V + tW + Z?') = -8A/do, 
f3" + {3\d ol ' + {3') = -8A/do, 

(d - l)i 2 + 1rf& = -8A/do. (22) 
Excluding /?', we get the equation for 7 

7 / ' + ^K + l)7 /2 + ^ = (23) 

Its solutions are different for positive and negative A. For A > we get (requiring 7(0) = 7'(0) = 
— a regular center at 1 = 0) 



e (do+i)7 = cos 2 (Al/); Al = / 2A ^±1. (24) 

V «0 

The last equation (122]) then gives for r = e 13 : 

r 2 = r 2 e - {d °- lh sin 2 (Ai/), r = const, (25) 

where, choosing ro = 1/Ai, we can satisfy the regularity condition r'(0) = 1. We thus have a 
configuration with a regular center but with a singularity e 7 — > and r — > 00 as Z — > 7r/(2Ai). 

For the potential (IT7|) . this case corresponds to £ > — 1. 

For A < 0, corresponding to e < —1, we have, instead of (|2"1]) and (125]) H, 



e 7 = cosh 2 /W' +1 )(A 2 /), A 2 = y-2A^t±. (26) 

sinh(A 2 Z) , v 

r r °[cosh(A 2 /)]W»- 1 )/( d o+D' 

and again, choosing r = 1/A 2 , we can satisfy the regularity condition r'(0) = 1. 

So the configuration with unbroken symmetry is completely regular and extends from the regular 
center I = to I = 00 where the warp factor e 27 and the radius r are infinitely growing functions. 

3.2. Solution with = 0* = const 7^ 

In this case, Eq. ( 120]) leads to a constant radius r = r* , and Eq. ({15"]) gives the relation 

2 /r 2 = -V{(f>*)/d , (28) 

whence it follows V(0*) = A/x 2 < 0. 
For 7 we get from ( [161) 



d 2 7 ' 2 = -2A, e do7 = e ± ^ =M ', (29) 



3 For d = 4 these formulae reduce to those found earlier in [23]. 
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Table 1: Classification of regular (do + 2)-dimensional solutions for arbitrary V(<p) by the types 
of asymptotic behavior at the largest or infinite I. The columns denoted r, 0, and 7 show their 
final values. Attraction or repulsion is meant with respect to the center. The symbol rj means a 
minimum of V(<p) . 



Notation 


/ range 


r 


<\> 


V(0) 


7 


Asymptotic type 


AO 


R+ 


00 


= 


V(0) < 


00 


AdS, attraction 


Al 


R+ 


00 


V 


V(rf) < 


00 


AdS, attraction 


A2 


R+ 


00 








00 


power-law, attraction 


BO 


R 


= r* 




K < 


±00 


horizon at one end 


Bl 


R + 




0*^0 


K < 


—00 


horizon, repulsion 


B2 


R + 


r* 


0* ^0 


K < 


00 


attracting tube 


C 


(0, Q 








V(0) > 


const 


second center 



and the coordinate range is I G R. In particular, for the potential (fTTj) we have from (1211) and ( 1T51) 

r = 1/0 - fl = const, (30) 
e = -2d (l - ft) ft - (1 - A 2 ) 2 < 0, (31) 

a = h 2 v(^) = -^y do ft(i - ft)}- (32) 

The configurations with <fi = const 7^ and r = const are regular but evidently cannot be inter- 
preted in terms of a brane world. They only provide the asymptotic behavior of the "tube" solutions 
presented below. 



4. Regular string solutions. Classification and map in the (T, e) plane 

The possible types of regular solutions to our field equations have been classified in Refs. [17,18]. 
Table 1 represents this classification for the string case d\ — 1. Compared to [17], this table omits 
the solutions existing only for d\ > 1, but additionally includes the solutions (labelled AO and B0) 
from Sec. 3. 

In what follows, we will deal with the potential fTTBl . For this potential, Fig. [1] presents the 
location of different regular string solutions in the plane of parameters (r, e) . The map shows 
solutions with the field having a constant sign. Those with alternate signs of <fi will be discussed 
in Sec. 6. There are no regular string solutions at e > 0. 

In Fig.[l], the black line e*(r) is the upper boundary of the area of class Al solutions with r — > 00 
as I — > 00. The points on the black line itself and in the whole area — 1 > e > e^(T) correspond to 
class B2 solutions. Fine-tuned solutions with two symmetric regular centers (see Fig. [7]) are located 
along the blue line. The class of fine-tuned solutions with a horizon (Bl) is represented by the red 
line. 

Let us briefly describe the classes of solutions presented in the map, postponing the derivation 
of some important details of the curves to the subsequent sections. 
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A1 



Figure 1: Location of solutions in the plane of parameters (I\ e). The black line e*(r) is the upper 
boundary of the area of class Al solutions with r — > oo as I — > oo. The black line itself and 
the whole area — 1 > e > £*(r) correspond to class B2 solutions with r(oo) = const. Fine-tuned 
solutions with two symmetric regular centers (class C) are located along the blue line £i(r). The 
class Bl of fine-tuned solutions with a horizon is represented by the red line £h(T). 

A: Configurations with infinite r 

From Table 1, one can notice a close similarity between the vacuum solutions AO (where = 0) and 
Al. In fact, the gravitational field in both cases is mainly governed by the (negative) cosmological 
constant. In the limiting case |e| ^> 1, e < the role of the symmetry breaking potential V(<f>) is 
negligible. Then, as follows from (|14p . 7' is always positive, the warp factor e 27 grows, and gravity 
is attractive towards / = . These solutions with r — > 00 at large I exist without any fine tuning. 

As |e| decreases, the potential V((f>) becomes more and more important, and at e approaching 
some £*(r) the derivative r' — > 0, so that class Al solutions pass over to asymptotically cylindrical 
fine-tuned class B2 solutions. 

The main features of class Al solutions at large I are: 

— the scalar field (f)(1) tends to a minimum of V(4>); 
- the quantities and (3'{l) tend to the same finite positive constant, so that e 7 ^ ~ r(l) 
grow exponentially. 

An example of an Al solution, found numerically, is presented in Fig.[2l 

As to the A2 regime, this is an exceptional solution corresponding to the condition V(0) = 0, 
hence e = — 1. In this case [18] a regular integral curve, starting at I = with = 0, finishes again 
with <p — > as / — > 00 . The large I behavior of r and the warp factor e 27 in the resulting regular 
solutions is 

r»i, e do7 ~ I. 
B: Asymptotically cylindrical configurations 

For these "tube" solutions, it is easy to verify that Eqs. (j2U]) - (p?T]) hold at large I and, in particular, 
/—>/* = const where < / 2 < 1 . 

The vacuum solution B0, with r = r* , actually interpolates between the cylindrical asymptotics 
of Bl ( e 7 — > 0, a double horizon) and B2 ( e 7 — > 00 at I — > 00, i.e., gravitational attraction towards 
the center 1 = 0). 
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Figure 2: Al solution with the parameters: T = 0.7, £ = —0.9. (a) — scalar field, (b, blue) - 
/?(/), (b, red) - , (c)-r(O- 



Eq. (13 ip does not contain x and allows finding the range of the input parameter e for which 
"tube" solutions are possible. The dependence £(/*) is shown in Fig.[3J By Eq. f!3T]) . 

ZcLq — 1 

For d = 4, e min = -16/7 = -2.2857.... 

It also follows from (1311) and Fig. [3] that in the range — 1 > e > e m \ n there are two branches of 
the inverse function f*(e) . In the range > e > — 1 there is only one branch. Other limiting values 
are expressed in terms of by fT29l - (!32l) . and 



7 



xny/dofiil-f*). (34) 




Figure 3: The dependence £(/*) fl3~T|) in the infinite cylinder case 

Type B2 configurations occupy a whole area in the (e, T) parameter plane, whereas Bl solutions 
require fine tuning and are located on the red line in the map. 



C: Configurations with two centers 

As was argued above, type C solutions can be symmetric and asymmetric with respect to reflections 
I — > l c — I . Symmetric solutions require one fine-tuning relation, which corresponds to particular 
curves in the (e, T) plane. The curve describing solutions with a constant sign of is presented 
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in Fig. [I] (the blue line). Other symmetric configurations will be discussed below. Asymmetric 
solutions can only appear at discrete points in the parameter plane, and we will not mention them 
any more. 

5. "Tube" solutions: location in the parameter plane 

The upper boundary e*(r) of type A solutions, found numerically point by point for do = 4, is 
presented in Fig. H] by the red line and the circles. Any point in the area e < e*(r), < T < oo 
corresponds to a type A solution with / monotonically growing from zero at the center to unity. 
The function £*(r) decreases from zero at T = to a minimum with e = e m { n according to Eq. fl33|) 
and then grows tending to — 1 as T — > oo. 

In the range > e > — 1 , the "tube" (fine-tuned) solutions exist only precisely on the line £*(T) , 
which comprises a border between A and C type solutions. 

In the range — 1 > e > e*(r), there are cylindrical solutions without fine tuning. This area is 
located between the zones of A and C type solutions. 



H 

190 




Figure 4: The upper boundary e*(r) of class Al solutions with infinite r. The points forming the 
red line are found numerically. The asymptotic dependences for r < 1 fj49l) and T ^> 1 (j4"0]) are 
shown by the blue and brown curves, respectively. 

In the limits of weak and strong gravitational fields (small and large T , respectively), a numerical 
analysis of the field equations is hindered, and we have derived the function £*(r) analytically. The 
blue line in Fig.H] corresponds to T 1 , the brown one to T 3> 1 . 

5.1. Strong gravity: r ^> 1 

As follows from numerical analysis, the scalar field in "tube" solutions is small in this case, and 
the potential can be expanded in a series: V(<f>) ~ Vo + ^V "cj) 2 , dV/dcj) ~ Vq4>. The problem is 
completely determined by the two constants Vo = V(0) and V Q " = (d 2 V '/ 'dtp 2 ) ^ = o . Let us introduce 
a new parameter A and a new function ip : 




A = 4x 2 Vo, ip — K 4>- 



(35) 
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For the Mexican hat potential (jTTI) . Vo = i] 2 (e + l)/4 and V " = — 1, so that 

A = T(e + l) (36) 

(recall that T := K 2 rf). Eqs. (|I5]) - (lI5j) take the form 

7" + 7W + /?) = -X/(2d )+^ 2 /d , (37) 

P" + P'(d i + P') = ^ 2 (l/d - e- 2/3 )-A/(2rf ), (38) 

V>" + VW + /?0 = ^(e- 2/3 -l) (39) 



and depend on only one dimensionless parameter A . The boundary conditions are 

7(0) = y(0) = V(0) = 0, ( e7| i=0 = 1, /?'(oo) = ^(oo) = 0. 

The scalar field equation ( |39|) is homogeneous with respect to ^ and looks linear; however, the 
system as a whole is nonlinear, and we have a nonlinear eigenvalue problem. The parameter do 
being fixed, there is only one dimensionless parameter A , whose ground state eigenvalue is expected 
to be of the order of unity. Then according to (I36I) for T ^> 1 , the parameter e is very close to — 1 , 
and fo in (I3T]) is ~ <C 1. From (1391) it follows that r = e^ — > 1 at large /. Nontrivial solutions 
exist for discrete values of A, one of which, corresponding to a monotonically growing , is found 
numerically: 

A = -7.433..., d = 4. 
The asymptotic dependence £*(T) for T ^> 1 

e = -l + A/r, T>1 (40) 
is presented in Fig. H] by the brown curve. 

5.2. Weak gravity: r«l 

The case T <C 1 is more complicated. A numerical computation shows (and it is verified analytically) 
that \e\ is exponentially small as r — > 0. From (13T|) we see that 1 — /q ~ —e/(2do) <C 1, and the 
limiting value of the circular radius (1301) . r* « ^j2do/\e is very large as compared with the "core" 
radius ~ 1 . The equations simplify differently in the two cases r<r» and r ^> 1 . The solutions 



should coincide in the intermediate region 1 <C r <C r* . 

It is convenient, for T 1, to rewrite the field equations in terms of r = e 13 : 

7" = - V(^-^)-^[s + d-/ 2 ) 2 ], (41) 

r" = W^ 7 - r - ^ - I/' 2 r - I [ e + (1 - / 2 ) 2 ] r, (42) 

/" = "/' ( d o7'-^)+^ -/(I-/ 2 )- (43) 



For r C r, we see that 7' ~ T < 1, and the term with r y' 2 in (|42p can be neglected. In the 
vicinity of the center, in the terms ~ V , we can set r = / and omit e . Then Eq. ( 1421) reduces to 

r" = -r^-^-(i-/ 2 ) 2 /, r«i, Z«r„ 
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Figure 5: Example of a type B2 solution for T = 2 and e = —1.1. (a) shows the scalar field f{l) 
(b) depicts and (3'{l) by the red and blue curves, respectively, and (c) displays r(7). 



where fo is the solution of fjl3l) with 7' = 0, r = I, and the boundary conditions /(0) = 0, 
/(oo) = 1. With r'(0) = 1, after integration we get 

-l 1 rl 



i + r 



-/ 2 lnZ + 2 / dlM' Q \nl-^- I dl(l-f ) 2 l . 
Jo ia o Jo 



The integrals quickly converge for / 3> 1 , and in the intermediate region 1 <C I <C r* we have 



„/2 



I +21 ( -lnr + 2J 2 - A ) . I I. I < r < r 5 . 



where the integrals J\ and J2 are found numerically: 

/*oo ^00 

Jl=/ ldl{l-flf = l, J 2 = ciz/o/ohl/ = 0.2. 



= 1, 

JO 



In the region r> 1, Eq. ( 1421) reduces to 
r" = -r ( - + r^-r ) , r<l, l<r. 
Taking into account that r' = at r = r* , we find 



-r' 2 = T (lnr* — lnr H — r 2 — r 2 

2 V 4rf o 4d 



r < 1, 1< r. 



In the intermediate region 1 r r* we have 

rl _ or I 'I r<1 l< r<n . 



,1 2d , 1 

2r - In In r 

2 e 2 



(44) 



(45) 



(46) 



(47) 



(48 



We have taken into account that r* = y2do7Rl • 

Comparing (JHj) with (HBl) . we find the fine-tuning relation between e and V for asymptotically 
cylindrical configurations in the weak gravity limit: 



-2d exp( -I-4J 2 -1- ^ 



. 1 1 

-0.33d exp — — + — 

r rf 



e _i.7 e - 1/r , rf = 4, T < 1. (49) 
This asymptotic dependence £*(r) for T 1 is presented in Fig. @] by the blue curve. 
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5.3. Solutions in the range — 1 > e > £*(T) 

In the range — 1 > e > e*(r), there exist type B2 solutions (r — > r* < oo as / — > oo) without any 
fine-tuning relation between the parameters e and V . The asymptotic values of the scalar field ( /* ) 
and the radius (r*) at large / do not depend on Y : 



2 = rfp-1- y/K-l) 2 +(g + l)(2rf -l) 

; * 2d - 1 

2 2d - 1 



r 



4 + ^o-l) 2 + (e + ip„-l) 



An example of such a solution is shown in Fig. [5] for T = 2 and e = — 1.1 . The scalar field /(/) is 
shown in Fig.[5h, and j3'(l) are red and blue curves in Fig. [5b, and r(Z) is displayed in Fig.[5b. 



5.4. Solutions with a horizon 

As to type Bl configurations with a horizon, in which 7(7) linearly decreases as / — > oo, their 
location, found numerically, is shown in Fig. [1] by the red line e = £fc(r), corresponding to a 
certain fine-tuning relation. An example of such a regular solution, with the parameters T = 2, 
e = —0.233846, is presented in Fig.[6j 




Figure 6: A type Bl solution with a horizon for T = 2 and e = —0.233846. Left: the scalar field <fi 
(red) and the metric function 7 (blue). Right: the circular radius r. 

The near-horizon metric has the asymptotic form 

ds 2 = C 2 e- 2hl V ^dx' M dx u - dl 2 - r 2 JQ 2 , h = 7 '(oo). (50) 

The substitution e~~ hl = p (converting I = 00 to a finite coordinate value, p = 0) brings the metric 
( l50l) to the form 

do 2 

ds 2 = C 2 p 2 ri iiu dx ii dx v - -f^ - r 2 JVl 2 , p -> 0. (51) 



Thus p = is a second-order Killing horizon in the 2-dimensional subspace parameterized by t 
and p; the extra-dimensional circular radius squared remains positive. It is of the same nature as, 
e.g., the extreme Reissner-Nordstrom black hole horizon, or the anti-de Sitter horizon in the second 
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Randall-Sundrum brane world model. A peculiarity of the present horizon is that the spatial part 
of the metric, which at large / takes the form p 2 (dx) 2 , is degenerate at p = 0. The volume of the 
do -dimensional spacetime vanishes as I — > oo. And it will remain degenerate even if we pass to 
Kruskal-like coordinates in the (t,p) subspace. But the D -dimensional curvature is finite there, 
indicating that a transition to negative values of p (where the old coordinate / no longer works) is 
meaningful. 

Thus solutions with strings (and/or monopoles) in extra dimensions may contain second-order 
horizons, and the degenerate nature of the spatial metric at the horizon does not lead to a curvature 
singularity; moreover, the metric may be continued in a Kruskal-like manner. However, the zero 
volume of the spatial section makes the density of any additional (test) matter infinite at p = 0. 
To consider these solutions as describing viable configurations, one needs to take into account the 
back-reaction of ordinary matter. It will evidently destroy such a configuration. 



6. Solutions with two regular centers: location in the parameter plane 

Symmetric type C solutions with two regular centers are located on the (e, T) plane in the region 
> e > — 1 to the right from the red fine-tuning curve £*(T) in Fig. H] or, which is the same, to 
the right of the black line in the full map, Fig.[TJ The solutions are fine-tuned, i.e., located along 
certain lines £at(T) in this region, where N is the number of half- waves, N — 1 is the number of 
knots (zeros) of the scalar field / = <p/i]. The point is that /, just as the radius r, is zero at both 
centers, but / can change its sign. So there are several families of regular solutions with different 
numbers iV of half- waves, each family corresponding to a line £at(A) in the parameter plane. The 
blue curve in Fig. Q] depicts £i(A). 

6.1. Solutions without knots of the scalar field 

In solutions where / has a constant sign, all three functions f(l) , r(Z) and ^(l) reach their extremum 
values at the equator I = l eq . Setting in the first integral ffTBl 

/'(/eq) = r'(l eq ) = 7 '(U = 0, (52) 
we find a relation between f(l eq ) ='■ / eq and r(l eq ) =: r eq : 



2 fl 

r ^ = l7i _ n _ f2 ^2- ( 53 ) 

J eq 



It is convenient to use (1531) together with (152]) as boundary conditions and perform numerical 
integration from the equator to one of the centers. Then the three conditions f(l c ) = 0, r(l c ) = 
and r'(l c ) = 1 determine the values of / eq and l c and a fine-tuning relation e = £i(r). 

An example of a configuration with two regular centers is presented in FigJT] for T = 2 ; the 
fine-tuned value of e is £i(2) = —0.3326..., it belongs to the line e = £i(T) . 

The (blue) curve e = £i(T) in Fig. [T] has been obtained numerically. For small and large values 
of T, this fine-tuning relation can be derived analytically. 



£i(r) for weak gravity, T 1 

This derivation repeats the one for Eq. (jUJ) . The main difference is that now we get from (I53|) the 
value of r(l) at the equator 

r eq = r(Jeq) = VVWl \e\ < 1 (54) 
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Figure 7: Example of a configuration with two symmetric regular centers without knots of the scalar 
field. T = 2, the fine-tuned value £i(2) = —0.3326.... It belongs to the line e = £i(r). Left: r(l) 
— red, f(l) — blue; right: 



and use it instead of r*. Substituting (J54"|) into fHTI) . we have instead of 

1 



J2 



2 

\s\ 



r I In — — 2 In r %-r 2 



r < 1, l«r. 



In the intermediate region 1 < r < r m it should coincide with (144"1) . The resulting relation is 
1 1 - Ji 



e = — 2 exp 



r d 

£i(r) for strong gravity, T ^> 1 



4J 9 -0.9 e" 1/r , 



T < 1. 



(55) 



Numerical integration shows that for T ^> 1 the scalar field / remains small in the whole interval 
between the centers, while ip — xr lf is °f the order of unity. Introducing, as before, A (}3"6"1) and 
taking into account that / < 1, we again find Eqs. (137|) — (I391) . It is convenient to integrate them 
from the equator to a center and to use boundary conditions in the form 



^(Zeq) = VK) = T'C^) = 0, 



cq- 



2 



2 

cq 



2^q-A" 



(56) 



The three parameters ip eq , Z c — Z e q and A should be determined from the conditions 

V>(k) = 0, r(/ c ) = 0, r'(Z c ) = 1. 

Numerical integration results in A = 1.9.... 

Thus, for type C solutions, the desired fine-tuning relation in the strong gravity limit is 



e = -1 + 1.9/r. 



(57) 
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6.2. Odd and even scalar fields 

In the above solutions, the scalar field f(l) without knots is an even function. 

Since /(£) may change its sign between the centers, there are two possibilities. If the number 
of knots of /(/) is even, then /(/) is an even function, reaching an extremum at the equator, and 
f'(hq) = 0. On the contrary, /(/) with an odd number of knots is an odd function: /(/ eq ) = 0, and 
f'(l) is then an even function having an extremum at I = l eq . 

Numerical integration of Eqs. (ll3l) - (fT5l) in the case of an even number of knots can be performed 
with the same boundary conditions (1561) as without knots. The results are displayed in Figs.lHHTOl 

Fig.El shows a few solutions for the scalar field /(/) with 2 knots and the corresponding functions 
r(i). 




Figure 8: Solutions with two symmetric regular centers with 2 knots of f(l) for T = 0.5, 0.75, 1, 
2, 3, 4 and 5 in the order of decreasing amplitude. The corresponding fine-tuned values of e are 
-0.51275, -0.62765, -0.7019, -0.8365, -0.888, -0.9146, and -0.93115. Left: f(l). Right: r(l); the 
additional dashed curve corresponds to the limit r — > oo . 

Fig.[9] shows the function 7(7) in the whole range of I and (b) in a close vicinity of the equator 
for visual clarity, to demonstrate a minimum at the equator. Recall that 7 enters into the equations 
only via 7' and 7", and so, without loss of generality, we have put j(l e q) = in Fig.[9] (b). The 
larger is T, the deeper is the local minimum of 7 at the equator. Altogether, the gravitational 
potential has three minima: one at the equator and two others near the regular centers. 

The fine-tuning relation e^T) for solutions with two knots of f(l) , found numerically, is shown 
in Fig.fTOl Each curve in Figs. [8] and [9] corresponds to a point on this curve. 

In solutions with an odd number of knots of f(l) , we have f(l e q) = 0, and from ([TBI) we obtain at 
the equator f' 2 = (e + l)/2. For numerical integration of Eqs. (ll3l) - (fT5l) in this case, it is convenient 
to use the following boundary conditions: 

P'(h q ) = /(*«,) = VGeq) = 0, /5(/ eq ) = /3 eq , f'(l eq ) = y/{e + l)/2. (58) 

Then the three conditions f(l c ) = 0, r(l c ) = and r'(l c ) = 1 determine the values of /3 eq and 
l c — l eq as well as the fine-tuning relation (e = £2(r) in the case of one knot). As an example, in 
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Figure 9: The functions 7(7) (a) in the whole range of / and (b) in a close vicinity of the equator 
for the solutions with the same set of parameters as in Fig. [HI 



T 
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I 



Figure 10: The fine-tuning relation e^T) for solutions with two knots of f(l) between the centers, 
displayed in Figs. [8] and [9j 



Fig.[TT]we present a numerical solution with one knot for T = 2, corresponding to e — —0.71. The 
left graph shows r(l) (red) and f(l) (blue). "/(I) is depicted on the right graph. The function 
is symmetric relative to the equator and has two minima close to the centers. 

One can note that if we include into consideration configurations with angular deficits and 
excesses at the centers, then the existence of solutions with two centers is not restricted to particular 
lines in the (e, T) plane. There will be a whole area of such solutions, bounded by e = min(e*(r), —1) 
from below, and by the line e = Sh(T) from above. Among the nonsymmetric solutions of this kind 
with several knots of the scalar field one can find those with multiple local maxima and minima of 
7(Z) . Their possible connection with matter trapping and the hierarchy problem is discussed below. 
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Figure 11: Example of a solutions with two symmetric regular centers and an odd scalar field with 
one knot. T = 2, the corresponding fine-tuned value S2(T) = —0.71. The left graph shows r(l) 
(red) and /(/) (blue); 7(/) is depicted on the right graph. 



7. Matter in the background of global string configurations 

In this section we discuss the problem of trapping of classical point particles and test scalar fields 
by the gravitational field of the global string configurations described above. 



7.1. Classical particles 

The motion of classical particles in the bulk may be equivalently described in terms of geodesies or 
using the Hamilton- Jacobi equation. We here use the second approach. 

The Hamilton- Jacobi equation for a point test particle of (primary) mass mo in space-time with 
the metric ffTOl) is 



e-27(0 



fdS\ 2 fdS 



- 'dsV i /as x2 



m = 

\dt J \dx J J \dl J r 2 (l) \d9 1 



The metric is homogeneous with respect to all coordinates except /, and the action can be written 
in the form 

S = Et-px + S l (l) + M9, (59) 



where E is the particle energy, p is the particle momentum along the coordinates x % , i — 1, do — 1 , 
6 is the angular coordinate in the extra dimensions (note that dfl 2 = d9 2 for d\ — 1) and M is its 
conjugate angular momentum. The remaining unknown function Si(l) satisfies the equation 




9 m M2 2 



where p 2 = E 2 — . Zeros of the square root determine the turning points of classical motion. 

Consider a particle with M = 0, i.e., moving in the bulk along the coordinate I (strictly to or 
from a brane if the latter is located at fixed /). Classical motion is allowed where the square root 
is real. The turning points It are determined by the equation 



P 



2 e -27(* t )_ m 2 = 0> 
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If there is a minimum of 7(7) at some I = Iq, a classical particle with p 2 = Trig e 27 ' ! °' cannot 
move along the / direction and is trapped precisely at this minimum of 7. Particles with slightly 
larger p 2 can move between two turning points in the vicinity of Iq . If it is a global minimum of 7 , 
particles with any p 2 > m 2 , e 2 "^ ) are trapped. 

It can also be verified that particles with the same value of p 2 but M 7^ (moving in the 8 
direction) have a still more narrow range of motion along I . 

In particular, near the equator of a configuration with two centers and three half-waves of <fi 
(see Fig.EJ), the turning points of finite classical motion exist for 7(Z eq ) < 7 < j m , where 7 m is the 
maximum of 7(/). Denoting m 2 q = mg e 27 ' lcq ' , we see that a classical particle is trapped near the 
equator if its energy is restricted by 

m 2 eCL <p 2 <m 2 eq e 2 ^-^. 



It moves along the Minkowski coordinates as a free particle of mass m, 



cq • 



7.2. Scalar fields 

Consider a test scalar field x with the Lagrangian L x such that 

2L X = d A X *d AX - m 2 oX *X (60) 

in the background of our string configurations with the metric ( II 01) . Here, the asterisk as a su- 
perscript denotes complex conjugation and mo is the initial field mass. The x fi e ld obeys the 
Klein-Gordon equation 

d A {V99 AB d B x) + yfgmlx = 0, (61) 

where g = \ det(gAB)\ — e Mo7+2/3 . Taking into account the symmetry of the problem, we can take 
a single mode of Xi assuming 

X {x A ) =X{l)e- ip ^ e me , (62) 

where p^ = (E,p) is the (do = 4) -momentum along the brane and n is an integer. Then X(l) 
obeys the equation 

X" + (d o y + p)X' + (p 2 e~ 27 - n 2 e~ 2/3 - m 2 )X = 0, (63) 

where p 2 = p^ = E 2 — p 2 is the effective mass squared, observed on the brane. 

As a trapping criterion for a mode X , it is reasonable to require finiteness of the x fi e ld energy 
E x per unit area of the brane, 

E x = J Ttix^dOdl = 2nj T*\xWgdl < 00, (64) 
where 

nix] = \ [ e" 27 (£ 2 + f) + X' 2 + (n 2 e- 2 ? + m 2 )X 2 } (65) 

is the temporal component of the x fi e ld stress-energy tensor. One can notice that the validity of 
( 1641 automatically guarantees finiteness of the norm J x fgx*xdld9 of the x fi e ld considered as a 
quantum-mechanical wave function. 
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The finiteness of E x in the background of different regular configurations with infinite extra 
dimensions described above depends on the behavior of solutions to (1631) at small and large I . 

We begin with considering the x ne ld behavior near a regular center 1 = 0, which is common 
to all classes of regular configurations. At small I , we have e 13 = r ~ Z and 7 — > . Hence Eq. (163 j) 
takes the approximate form 

/X" + X' + Z(p 2 - ml)X = (n = 0), (66) 
IX" + X' - {n 2 /l)X = (n ^ 0). (67) 

Eq. (166]) is solved by zero-order cylindrical functions if p 2 7^ wig and in elementary functions if 
p 2 = m 2 ; Eq. ( 1671) is an Euler equation. The solutions behave as follows at small Z : 

X t* C 1 + C 2 hil (n = 0), 

X^C 3 r + C 4 r n (n^O), (68) 

with integration constants Cj. To make the integral in (1(31) converge as Z — > 0, one should choose 
C2 = and C4 = 0, i.e., in each case, only one of the two linearly independent solutions. 

Now, consider the asymptotic form of solutions to Eq. (163j) as I —>■ 00 for different background 
configurations. 

Al: at large I, 7 <~ (3 ~ hi, h = const > 0. In Eq. (1631) . the terms with p 2 and m? are negligible, 
and the solution has the asymptotic form 



X ^C + e~ a+l + C_e- a ~ l , 2a± = (D-l)h±^(D-l) 2 h 2 + 4m 2 , (69) 

where C± are integration constants and D = do + 2 is the full space-time dimension. It is easy to 
verify that the criterion (164|) holds for the solution with C + 7^ 0, C_ = 0. Thus scalar fields with 
any nonzero mass can be trapped on such branes. 

A2: at large I, e d ° 7 ~ e 13 ~ I. Again, the terms with p 2 and m 2 are negligible, and Eq. ( |63l 
transforms to X" + IX' j I — m^X = whose solution is 

X w (C+ e moi + C_ e- mo/ )/^, (70) 
and evidently the solution with C + = satisfies the criterion f l64l) . 
Bl: at large /, r = e' 3 — > r*, 7 ~ — /iZ, h > 0, and the approximate form of Eq. fl63|) is 

A // -rf / i A / +p 2 e 2W A = 0. (71) 
For p 7^ , it is solved in cylindrical functions, the general solution being 

X = e W%7 2 (Me M ) ~ e ^- 1 ) h '/ 2 sin(M e ^ + $), (72) 

where $ is a constant phase. It is easy to verify that E x diverges as j e hl dl . So massive modes 
with any p 2 > are not trapped by Bl configurations. 

B2: at large I, r — * r* and 7 w /iZ, Zi > 0. The situation is almost the same as in case Al; the 
solution to (16"31) has the asymptotic form (1591) with the replacements 

D — 1 1— > d , wiq + n 2 /r 2 . 

Again, only the solution with C_ = provides convergence of E x . 

Thus the configurations of classes Al, A2 and B2 can trap massive scalar modes; since at both 
large and small Z only one of the two linearly independent solutions to (|63|) is selected, we have 
a boundary- value problem with a discrete spectrum of p 2 for any given values of mo , n and the 
background parameters. 
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7.3. The Schrodinger equation 



It is helpful to reformulate the boundary-value problem for scalar field modes in terms of the 
Schrodinger equation. To do so, we make the following substitutions in (|63|) : 



dl = e^dx, 



X(l) = y(s)/v7to 



(73) 



where f(x) = e^ " 1 ^ 3 . The new variable x is actually an analogue of the well-known "tortoise 
coordinate" in the analysis of spherically symmetric metrics, such that the metric takes the form 



ds 2 = e^^dx^dx" - dx 2 ) - e 2p d6 2 . 
Then Eq. (1631) transforms to the Schrodinger form 

Vxx + [p 2 - V eS (x)]y = 
with the effective potential 



cff 



/ 2 i 2 -20\ 27 , fxx fx 

(m + n e 1 ) e 1 + -^j 



2f 4/2' 



(74) 



(75) 



(76) 



where the subscript x denotes d/dx. Recall that the eigenvalue p 2 is the effective mass squared, 
observed in Minkowski space. 

Near the center (without loss of generality, x ~ I — > 0), we have 



cff 



n 2 /x 2 + m 2 + - [1 + 2(do - l)^ xx + (3 X 



J x=0 



(77) 



So for n 7^ it is a potential wall whereas V e g — > const for n = . 
At large I for different backgrounds, we have: 



Al : 



A2 



Bl 



B2 



x 

V eS (x 



x + < oo; x + — x ~ e 
h 2 



-hi 



h>0: 



,2W 



m 2 + —(d 2 + 2d ) 



oo; 



X ~ l^-V/d* 

V eS {x) « m 2 e 2 ^~ / 2 / rf0 ~ x 2 ^- 1 ). 
e M -> oo, /i > 0; 



Kfr(x) 



-2M 



n 



2 h 2 



(78) 
(79) 
(80) 



x 



< oo; x* — x ~ e 



/i > 0; 



2hl 



n 2 h 2 



We see that in the cases Al, A2, B2 the potential is rising to infinity at large I, which leads 
to discrete spectra of p 2 . For Bl configurations (with a horizon at I = oo), in usual quantum 
mechanics we would expect a continuous spectrum of states; in our case, with the appropriate 
boundary conditions, as we saw above, there are no admissible states with p 2 > . 
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7.4. Massless modes in configurations with infinite extra dimensions 

For a possible massless mode, p 2 = ml = n 2 = 0, Eq. (1631) is easily solved: 

I^de-* 7 ^, X = Cl J e ~ d0l ~ Pdl + C * Ci, 2 = const. (82) 

and X is found by quadrature. 

One of the solutions is X = const. It is easy to verify that with this solution, which is well- 
behaved at a regular center, the energy E x diverges at large I in the backgrounds Al, A2, B2 
but converges in the background Bl. 

For the other solution with C\ ^ , on the contrary, E x converges at large / in the backgrounds 
Al, A2 and B2 and diverges in Bl. This solution, however, is singular at the center and leads there 
to a divergence in E x . 

Thus Bl configurations with horizons, being unable to trap massive scalar fields, are the only 
ones that can trap a massless scalar. 



8. Configurations with two centers and the hierarchy problem 

In configurations with two symmetric regular centers and two knots of the scalar field, there are 
three minima of the "gravitational potential" 7, see Fig.|9j The minimum at the equator is higher 
than the other two, located near the centers. A similar (and even more complicated) structure may 
be expected for configurations with a larger number of scalar field knots. The minima of 7 are able 
to trap classical particles. As to quantum ones (at least for zero spin), the effective potential (175]) 
not necessarily has a minimum precisely where 7 has a minimum, and an additional detailed study 
is necessary. Though, at least quasiclassically, quantum and classical particles must be trapped in 
close positions, and the main difference between them is that quantum particles can tunnel from a 
higher minimum of V e R to a lower one. 

Now, suppose a particle described by a certain mode of the \ ne ld (for simplicity, with n — 0) 
is trapped at some position U . The mode equation (|63|) may be rewritten in the form 

(y/gX')' + e"V = ^m 2 X. (83) 

Let us integrate this equation over the extra dimension from one center to the other. The term 
/ [yfgX'j dl = because yjg = re* 7 is zero at both centers. For a particle trapped at some fixed 
position I = li, we get 

To interpret this result, we note that the whole picture looks quite different depending on the 
size of the extra dimensions, characterized by the distance l c between the centers. This size, in turn, 
varies with the value of T = x 2 7/ 2 : it is close to unity (i.e., the length unit which is also arbitrary) 
for large T and tends to infinity as r — >• . For (comparatively) weak gravity of the string, r — > , 
when l c is very large, all minima of ^(l) form individual branes located in the bulk far from one 
another. In this case, an observer located on one of the branes sees only particles corresponding to 
modes trapped on this brane; tunnelling from one brane to another will be seen as appearance or 
disappearance of observable particles. The whole picture may be used for treating the interaction 
hierarchy problem in the spirit of Randall and Sundrum's first model [4] . 
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In the opposite case r — > oo (and provided the unit length (Ao?7)~ 1//2 is also sufficiently small), 
l c can be a length invisible for modern instruments, say, l c <C I0~ 17 cm. We then arrive at a 
picture close to the original Kaluza-Klein concept; particles with the same primary mass m , being 
trapped at different minima of the effective potential, are seen as particles with different masses, 
and the tunneling process from a higher minimum to a lower one is observed as a decay of a particle 
of larger mass to that of smaller mass, with energy release in some form. This may be a natural 
explanation of the existing families of particles with different masses but similar other properties. 
A more detailed study of this opportunity is desirable but is beyond the scope of the present paper. 

9. Conclusion 

Our phenomenological approach based on the macroscopic theory of phase transitions with sponta- 
neous symmetry breaking allows us to study the general physical properties of topological defects in 
the framework of the brane world concept. In particular, in this paper we have studied the gravita- 
tional properties of global strings located in extra dimensions. We have given a general description 
and classification of possible regular solutions and presented a map showing the location of different 
solutions in the space of physical parameters. 

Among the variety of regular solutions, there are ones possessing brane features, including those 
with multiple branes, as well as those of potential interest from the viewpoint of the hierarchy 
problem. 

In connection with branes, we have analyzed the possibilities of trapping of classical particles and 
scalar fields. We have shown that, contrary to a domain wall, in the case of an extra-dimensional 
global string, matter can be trapped by gravity even without coupling to the scalar field that forms 
the string itself. 

Among the configurations with two centers, the structures having several minima of 7(7) may 
be of interest in connection with the hierarchy problem. If the distance between the centers is small, 
we work within the Kaluza-Klein concept, and the same particle, being trapped at different minima, 
looks for an observer as a family of similar particles with different rest masses. 
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